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Abstract 

We evaluate all two-body decay modes of the gluino, in the Minimal Supersymmetric 
Standard Model with complex parameters (cMSSM). This constitutes an important 
step in the cascade decays of SUSY particles at the LHC. The evaluation is based on a 
full one- loop calculation of all two-body decay channels, also including hard QED and 
QCD radiation. The dependence of the gluino decay to a scalar quark and a quark on 
the relevant cMSSM parameters is analyzed numerically. We find sizable contributions 
to the decay widths and branching ratios. They are, roughly of 0(zb5%), but can go 
up to ±10% or higher, where the pure SUSY QCD contributions alone can give an 
insufficient approximation to the full one- loop result. Therefore the full corrections are 
important for the correct interpretation of gluino decays at the LHC. The results will 
be implemented into the Fortran code FeynHiggs. 
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1 Introduction 



One of the most important tasks at the LHC is to search for physics effects beyond the 
Standard Model (SM), where the Minimal Supersymmetric Standard Model (MSSM) [1] 
is one of the leading candidates. Supersymmetry (SUSY) predicts two scalar partners for 
all SM fermions as well as fermionic partners to all SM bosons. Especially it predicts the 
fermionic gluino as a superpartner of the bosonic gluon. 

If SUSY is realized in nature and the gluino and/or scalar quarks are in the kinematic 
reach of the LHC, it is expected that these strongly interacting particles are copiously pro- 
duced Refs. [2],[3]. After the (pair) production of these particles they are expected to decay 
via cascades to lighter SUSY particles and quarks and/or leptons. Many models predict 
that the heaviest SUSY particle is indeed the gluino (see, for instance, Ref. |1] and refer- 
ences therein), and consequently, the decay to a scalar quark and a quark will one of the 
most relevant decays in such a cascade. Only if these two-body decays are kinematically 
forbidden other decay channels can become numerically relevant. These are either two-body 
decays that are purely loop-induced, or three-body decays via a virtual squark. In this paper 
we will concentrate on the two-body decays that are allowed at tree-level already. We will 
also, for the sake of completeness, evaluate the two-body decays that appear only at the 
loop-level, the decays of the gluino to a gluon and a neutralino, which can potentially serve 
as a production source of the cMSSM cold dark matter (CDM) candidate. 

In order to yield a sufficient accuracy one-loop corrections to the various gluino decay 
modes have to be considered. As outlined above, we take into account all two-body decay 
modes of the gluino in the MSSM with complex parameters (cMSSM). More specifically we 
calculate the full one-loop corrections to 



A difference between r(^ — )■ qj q) and r(^ — )■ g) can arise from complex parameters. In 
the case when all parameters are real (and sometimes to simplify the notation) we also use 
the notation r(^ — )■ q) referring to both channels. The total width is defined as the sum of 
the channels ([I]) - ©j where for a given parameter point several of the 28 channels may be 
kinematically forbidden. As it is expected, the channel (E]) is numerically irrelevant for the 
total width as long as one or more of the two-body tree-level decays is kinematically allowed. 

One-loop QCD corrections to gluino decays have been evaluated in Ref. [5]. The calcu- 
lation was done in the MSSM with real parameters (rMSSM). Those 0{as) corrections have 
been implemented into the code SDECAY [6] . To our knowledge, no evaluation of electroweak 
(EW) corrections to gluino decays, nor any higher-order corrections involving complex phases 
have been performed so far. Gluino decays at the tree-level were employed, for instance, in 
Ref. [7] to determine their Majorana/Dirac character. Gluino polarizations were analyzed 
in Ref. [8]. Again, both analyses were performed in the rMSSM. 

Several methods have been discussed in the literature to extract the complex parameters 
of the model from experimental measurements. However, no such an analysis, to our knowl- 
edge, of the phase of the gluino mass parameter, M3, has been performed so far. We analyze 



r(^ ^ qjq) 



(z = 1, 2; q = t,b,c, s, u, d) , 
(2 = 1,2; q = t, b, c, s, u, d) , 
(A; = 1,2,3,4) . 



(1) 

(2) 
(3) 
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the effects of tliis pliase on tlie brancliing ratios (BR's) of tlie gluino, also as a motivation 
to devise experimental strategies to determine this parameter at the LHC (or other future 
colliders) [H]. 

In this paper we present for the first time a full one-loop calculation, including electroweak 
effects, for all two-body decay channels of the gluino in the cMSSM. The calculation includes 
soft and hard QED and QCD radiation. In Sect. [2] we review the relevant sectors of the 
cMSSM. Details about the calculation can be found in Sect. [3l and the numerical results for 
all decay channels are presented in Sect. HI were we discuss especially the size of the EW 
corrections and the effects from the complex gluino phase. The conclusions can be found in 
Sect. [5l The evaluation of the branching ratios of the gluino will be implemented into the 
Fortran code FeynHiggs [T0HT3] . 



2 The relevant sectors of the complex MSSM 

All the channels ([1]) - (jS]) are calculated at the one-loop level, including hard QED and QCD 
radiation (for the decays that exist already at the tree-level). This requires the simultaneous 
renormalization of several sectors of the cMSSM, where details can be found in Refs. [TIHTB] . 
In the following subsections we briefly review these relevant sectors. 



2.1 The squark sector of the cMSSM 

For the evaluation of the one-loop contributions to the decay of the ^ to a squark and quark, 
a renormalization of the scalar quark sector is needed. The squark mass matrix read 

^/Ml+ml + Ml ( J3 - Q^sl) cos 2/3 m,X* \ 

' V ^9^9 Ml + ml + MlQ,slcos2pJ 

with 

Xq = Aq — jj* n , K = {cot /3, tan/3} for q = {m— type, rf— type} . (5) 

Mq^ and Mq^ are the soft SUSY-breaking mass parameters, where M^^ is equal for all 
members of an SU{2)l doublet and Mq^ depends on the scalar quark flavor, rriq is the 
mass of the corresponding quark. Qq and Iq denote the charge and isospin of g, and Aq is 
the trilinear soft-breaking parameter. Mz and My/ are the masses of the Z and W boson, 
Cw = Mw/Mz, and s„ = ^yl — c^. The mass matrix can be diagonalized with the help of a 
unitary transformation Ug, 

The mass eigenvalues depend only on \Xq\. The scalar quark masses will always be mass 
ordered, i.e. rriq-^ < rriq^: 

= \{Ml + Ml) +ml + \ll Ml cos 2/3 
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T -Ml + Ml {II - 2Q,sl) COS2/3]' + Aml\X,\^ . (7) 

Details about the renormalization of the scalar quark sector of the cMSSM can be found in 
Refs. p^[T6] . The employed renormalization preserves the SU{2)l symmetry of the model. 
It is ensured that all external particles fulfill the required on-shell (OS) properties for masses 
and Z factors. The field renormalization requires the renormalization of the off-diagonal 
entry in the squark mass matrix, which leads to a renormalization of the trilinear couplings 
Aq, /i and tan/3, see below. It was shown in Refs. [T5l[T6] that the renormalization produces 
stable and well-behaved results for nearly the whole cMSSM parameter space. 

Finally it should be noted that we take into account the absorptive part of the self-energy 
type contributions on the external legs via combined Z factors which are different for incom- 
ing squarks/out going antisquarks (unbarred) and outgoing squarks/incoming antisquarks 
(barred) [16]: 

(i) The diagonal Z factors read 

l^2^u=\^2^u = -^.HJ (^ = 1,2). (8) 

(ii) The off-diagonal Z factors read 

«Ji2 (^i^_^2j "J 21 ^ml-ml) 

Hg-. denotes the scalar quark self-energies, and H'^m^) = 1^2-^2 • See Ref. [16] for the 
other renormalization constants and further details. 



2.2 The quark sector of the cMSSM 

In this section we briefly describe the quark sector of the cMSSM and its renormalization, 
extending the corresponding discussions in Refs. [T514T7] . The quark mass, m^, and the quark 
fields ql, Qr are renormalized in the following way: 

rriq mq + 6mq , (11) 
Ql/r ^ (1 + l^^^/"") Ql/r , (12) 
qL/R^il + lSZ^^^)qL/n, (13) 

with 6mq being the quark mass counterterm and 6Zq being the combined Z factors of the 
left /right-handed quark fields, respectively. They are determined separately to include effects 
from the absorptive parts of the self-energy type contribution on an external quark leg [16] . 
The unbarred quantities denote incoming quarks/outgoing antiquarks and the barred denote 
outgoing quarks/incoming antiquarks. The renormalized self energy. Eg, can be decomposed 
into left /right-handed and scalar left /right-handed parts, E^^^ and Tig^^^^, respectively, 

tq{p) =^co^t^ip')+ ^u+tfip') + a;-Sf (p^) + co+t'/ip') , (14) 
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where the components are given by 



YSR(J2\ 



(5m, 



q 1 



dm 



Q 1 



(15) 
(16) 
(17) 



and cj± = i(l ± 75) are the right- and left-handed projectors, respectively. 

The quark mass is defined on-shell [18], yielding the one-loop counterterm ^m^: 



5m 



OS 



|Re [niq [S 



J(mJ) + (mj 



f K) + (m; 



} 



:i^ 



referring to the Lorentz decomposition of the self energy Sg(p), see Eq. (IT^ . Re denotes 
the real part with respect to contributions from the loop integral, but leaves the complex 
couplings unaffected. 

Special care is needed for the bottom sector due to potentially large effects at large 
tan/3 (denoting the ratio of the two vacuum expectation values). If (and only if) there are 
no external bottom quarks, the bottom-quark mass is defined DR (see Ref. [16] for more 
details), yielding the one- loop counterterm 5m}^^: 



iRe {mb [S^ 



SL, 



SRi 



(19) 



The new (diagonal) field renormalization constants, taking into account the absorptive 
part of the self-energy type contribution on an external quark leg are different for incoming 
quarks/outgoing antiquarks (unbarred) and outgoing quarks/incoming antiquarks (barred). 



± 



1 



2mr. 



fK) 



1 



2 m, 



f(-;) 



sf(m;) 



Ef(m^ 



)] 



(m^) + Sf 



- S^^fm? 



m„ 



m„ 



+ ma 



Sf'(mJ) 



SR' 



■^SR' 



m„ 



m„ 



(20) 



(21) 



For further details see Ref. [T6]. 



The input parameters in the b sector have to correspond to the chosen renormalization. 
We start by defining the bottom mass, where the experimental input is the SM MS mass [T9] . 
The value of m^^(/i/j) (at the renormalization scale fin) is calculated from m^^(mb) at the 
three loop level following the prescription given in Ref. [20]. The input parameters can be 
found in Sect. 14.11 The "on-shell" mass is connected to the MS mass via 



m 



MS 




(22) 
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The DR bottom quark mass is calculated iteratively from 





1 + Aft 


1 + 6m^^ - 5mP^ 




1 + Ab 





(23) 



with an accuracy of |1 — (m^^)'^"V('^6*^)''""^''l < reached in the nth step of the iteration 
and Sm^^ as given in Eq. (fTSl) with q = b. The quantity A;, [2T1I22] resums the 0{{as tan/3)") 
and O {{at tan (3)"') terms and is given in Eq. (66) of Ref. [T6] . 

In the case of external bottom quarks we use an OS renormalization scheme for the 
bottom sector. The bottom quark mass is then obtained from 

mf^=mf^ + 6mf^-6m^\ (24) 

again with Sm^^ from Eq. f|T8|) with q = b. 



2.3 The gluino and the strong coupUng constant 

The soft-breaking gluino mass parameter M3 is in general complex, 

M3 = iMsle*'^^ (with the gluino mass m-g = iM-s]) . (25) 

We choose an OS renormalization for the gluino, where the renormalization constants can 
be found in Ref. [H]. It properly takes into account the complex phase of M3 and yields 
an on-shell gluino, as it is required for an external particle in our decaysQ It is possible to 
modify the field renormalization constants to take into account the absorptive part of the 
self-energy type contribution on the external gluino leg. The renormalization constants then 
read (TB] (unbarred (barred) for an incoming (outgoing) gluino), 

SM, = % {m, [Sf (m|) + Sf (m|)] + [sf (m^) + Sf (m|)] } e^^^ , (26) 



^^(-i) + + (m^)) + m, (Sf (mi) + Sf (^1; 

6Z^^^ = 5Zf'^ . (28) 

The last formula holds due to the Majorana character of the gluino and we have chosen 
5ipg = 0. This choice is possible, since the imaginary part of M3 does not contain any 
divergence 



Our renormalization of the strong coupling constant, as, is described in Refs. 
(and references therein). The decoupling of the heavy particles is taken into account in the 
definition of a^: Starting point is a^^{Mz) [I9], where the running can also be found in 



^ The general renormalization procedure is described in Refs. [T4lll6j . 

^ Formulas including the Sipg contributions can be found in Sect. 2.1.2 in Ref. 
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Ref. [19]. From the MS value the DR value is obtained at the two-loop level (with rij = 6 
for jipi> rrit) via 



R) 



MS,(6)/ N ^7 / MS,(6)/ , N\2- 



Within the MSSM as to one-loop reads 



47r 



DR,(6)/ N / 

1 + ^^ ^fln^ + ln^ 



8 TT^ 



(29) 



(30) 



with Mq = Ilq{mq-^mq^) . The log terms origin from the decoupling of the SUSY QCD 
(SQCD) particles from the running of ag at lower scales fiR < /idee. = "^g- 



2.4 The Higgs and gauge boson sector of the cMSSM 

The two Higgs doublets of the cMSSM are decomposed in the following way, 



Hi 

n2 



H22 



-0r 



V2 + 7|(02 + 



(31) 



Besides the vacuum expectation values Vi and V2, in Eq. (13T]) a possible new phase ^ between 
the two Higgs doublets is introduced. The Higgs potential Vh can be written in powers of 
the Higgs field^, 



V, 



H 



2 (01,02,Xl,X2) M^^^^ 



<P2 

Xi 

\X2j 



(32) 



where the coefficients of the linear terms are called tadpoles and those of the bilinear terms 
are the mass matrices M,^^^^ and M^±^±. After a rotation to the physical fields one obtains 



V, 



H 



Thh + TuH + TaA 



-{h,H,A,G) Mt-\^ 



H 

A 

\GJ 



(i/+, G+) 



G- 



+ ... 



(33) 



Corresponding to the convention used in FeynArts/FormCalc, we exchanged in the charged part the 
positive Higgs fields with the negative ones, which is in contrast to Ref. [13]. As we keep the definition of 
the matrix M^i^i used in Ref. [13] the transposed matrix wih appear in the expression for M^^^g,^. 
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where the tree-level masses are denoted as m/j, rriH, tra, rnc, Mh±, t^g^- With the help 
of a Peccei-Quinn transformation [21] fi and the complex soft SUSY-breaking parameters in 
the Higgs sector can be redefined [23] such that the complex phases vanish at tree-level. As 
input parameter we choose the mass of the charged Higgs boson, Mh±. All details can be 
found in Refs. 

Higgs bosons, hP , H^, H^, , and the electroweak gauge bosons, Z, W^, 7, appear 
only as internal particles, hence no renormalization is required. Furthermore we use tree- 
level masses and couplings for the Higgs bosons that appear internally in the loops, see 
Sect. [3l The only exception is the renormalization of tan (5 that enters the squark field 
renormalization. Here we choose a DR scheme as defined, for instance, in Eq. (120c) in 
Ref. m. 



2.5 The chargino/neutralino sector of the cMSSM 

The mass eigenstates of the charginos can be determined from the matrix 

VV2cos/3MvK At / ■ 

In addition to the higgsino mass parameter it contains the soft breaking term M2, which 
can also be complex in the cMSSM. The rotation to the chargino mass eigenstates is done by 
transforming the original wino and higgsino fields with the help of two unitary 2x2 matrices 
U and V, 

X- = with # = t/,, i^l^ and = V,, i^l^ , (35) 

where the mass eigenstate can be expressed in terms of either the Weyl spinors ipf and 
-0/^ or the Dirac spinor Xi with i, j = 1, 2. 

These rotations lead to the diagonal mass matrix 

M,-=V*X^Ut= I \ . (36) 

From this relation, it becomes clear that the chargino masses m-± can be determined as the 
(real and positive) singular values of X. The singular value decomposition of X also yields 
results for U and V. 

A similar procedure is used for the determination of the neutralino masses and mixing 
matrix, which can both be calculated from the mass matrix 



Y 



/Ml -Mzs^cos(3 Mzs„sm/3 \ 

M2 Mzc^ cos (3 -MzCw sin /3 

-Mz cos (3 Mz Cw cos (3 — /i 

y Mz sin /3 -Mz sin /3 -/i J 



(37) 
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This symmetric matrix contains the additional complex soft-breaking parameter Mi. The 
diagonalization of the matrix is achieved by a transformation starting from the original 
bino/wino/higgsino basis, 



xl 



]¥"> 



N*YN 



/m^o 
m^o 

A2 





V 











A3 





\ 







(3^ 



where ip^ denotes the two component Weyl spinor and Xfc the four component Majorana 
spinor of the kth neutralino field with k,l = 1,2,3,4. The unitary 4x4 matrix N and the 
physical neutralino masses again result from a numerical singular value decomposition of Y. 
The symmetry of Y permits the non-trivial condition of using only one matrix N for its 
diagonalization, in contrast to the chargino case shown above. 

A renormalization of this sector, as described in detail in Refs. [IU[T6], is not required. 
Neutralinos and charginos appear only as internal particles in the loop corrections of the 
two-body decays, see Sect. [31 In addition, neutralinos appear as external particles only in 
purely loop- induced processes, see Eqs. ([I]) - ([3]). Consequently, we obtained our results by 
using tree-level masses and couplings in the chargino and neutralino sector. Again, the only 
exception is the renormalization of fi entering the squark field renormalization. We follow 
the on-shell prescription given in Ref. [I6], where 6fi is defined in Eq. (181). 



3 Calculation of loop diagrams 

In this section we give some details about the calculation of the higher-order corrections to 
the gluino decays. Sample diagrams are shown in Figs. [H O Not shown are the diagrams for 
real (hard and soft) photon and gluon radiation. They are included via analytical formulas 
following the description given in Ref. [18] . 

The internal generically depicted particles in Figs. [11 [2l are labeled as follows: F can be 
a quark, chargino, neutralino or g, S can be a squark or a Higgs boson, V can be a 7, Z, 
or g. Internally appearing Higgs bosons do not receive higher-order corrections in their 
masses or couplings, which would correspond to effects beyond one- loopfl 

In Fig. [Hwe have furthermore omitted in general diagrams of self-energy type of external 
(on-shell) particles. While the real part of such a loop does not contribute to the decay width 
due to the on-shell renormalization, the imaginary part, in product with an imaginary part 
of a complex coupling (in our case coming from (fg) can give a real contribution to the decay 
width. While these diagrams are not shown explicitly, they have been taken into account 
in the analytical and numerical evaluation via renormalization constants, see Sect. [2l The 
impact of those absorptive contributions will be discussed in Sect. 14.31 

The diagrams and corresponding amplitudes have been obtained with FeynArts \26\- The 
model file, including the MSSM counter terms, is based on Ref. [HHTB] . i.e. to match exactly 
the renormalization prescription described in those articles (see also Ref. [27] for other sec- 
tors). The further evaluation has been performed with FormCalc (and LoopTools) [28]. As 

We found that using loop corrected Higgs boson masses for the internal Higgs bosons also result in 
finite results, where the numerical difference for the various decay widths is negligible. 
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Figure 1: Generic Feynman diagrams for the decay g ^ qtq {i = 1,2). F can be a quark, 
chargino, neutralino or g, S can be a squark or a Higgs boson, V can be a 7, Z, or g. 
The tree-level diagram is not shown here. 




Figure 2: Generic Feynman diagrams for the decay (7 — )■ = 1? 2, 3, 4). Here F denotes 

a quark and S a squark. It should be noted that this process exists only at the loop-level. 



regularization scheme for the UV-divergences we have used constrained differential renor- 
malization [29], which has been shown to be equivalent to dimensional reduction [30] at 
the one-loop level [28]. Thus the employed regularization preserves SUSY [311|32]. All UV- 
divergences cancel in the final result. 

The IR-divergences from diagrams with an internal photon or gluon have to cancel with 
the ones from the corresponding real soft radiation. In the case of QED we have included the 
soft and hard photon contribution following the description given in Ref. [18]. In the case of 
QCD we have modified this prescription by replacing the product of electric charges by the 
appropriate combination of color charges (linear combination of Ca and Cp times a^). The 
IR-divergences arising from the diagrams involving a 7 or a (7 are regularized by introducing 
a finite photon mass, or gluon mass, A. While for the QED part this procedure always 
works, in the QCD part due to its non-abelian character this method can fail. However, 
since no triple or quartic gluon vertices appear, A can indeed be used as a regulator. All 
IR-divergences, i.e. all divergences in the limit A — ?■ 0, cancel once virtual and real radiation 
for one decay channel are addedl 

^ Using tree- level, as well as one- loop down- type squark masses (see Sect. 12. ip . yields a cancellation of 
IR divergences to all orders for all g decays. 
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For completeness we show here also the formulas for the tree-level decay widths: 



2^ 



+ 4 Re{C(^, g„ q)l C{g, g„ g)^} m, ^ ' (39) 

where \{x,y,z) = [x — y — z)"^ — Ayz and the couplings C(a, 6, c) can be found in the 
FeynArts model files [33]. C(a, 6, 0)2,^^2 denote the part of the coupling which is proportional 
toi(llT75). 

4 Numerical analysis 

In this section we present a numerical analysis of all 28 decay channels ("xy"). In the various 
figures we show the decay width and its relative correction at the tree-level ("tree") and at 
the one-loop level ("full"), 

"pfuU "ptree 

The total decay width is defined as the sum of the 24 decay width^ , 

■pfuU ptree 

Tr::^Y.^'-{g^^y) . rf"i^ = ^rf"^^->xy), = ^^^^^ • ^^^^ 

xy xy *°t 

We also show the absolute and relative changes of the branching ratios, 



r^""(^^xy) RRfuu _ r^-n^ ^ xy) _ biT^^br^ 

J- tot ^ tot Dti. 



The last quantity is crucial to analyze the impact of the one-loop corrections on the phe- 
nomenology at the LHC. 

4.1 Parameter settings 

The renormalization scale, /i/j, has been set to the mass of the decaying particle, i.e. /ir = rrig. 
The SM parameters are chosen as follows, see also 

• Fermion masses: 



rUe 


= 0.51099891 MeV , 




= MeV , 




= 105.658367 MeV , 




= MeV , 




= 1776.82 MeV , 




= MeV , 




= 53.8 MeV , 


rrid 


= 53.8 MeV , 



^ We neglected the g ^ xlg decays in the total decay width due to their negligible contribution, see 
Sect. IMI and Sect. 1431 
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rUc = 1.27 GeV , 
mt = 172.0 GeV , 



= 101 MeV , 
mp(mb) = 4.25 GeV . 



(43) 



THu and are effective parameters, calculated through the hadronic contributions to: 

f=u,c,d,s,b J 



The CKM matrix has been set to unity. 
Gauge boson masses: 

Mz = 91.1876 GeV , Mw = 80.399 GeV , (45) 



Coupling constants: 

a = 1/137.035999679 , af^{Mz) = 0.1184 . (46) 



The Higgs sector quantities (masses, etc. at the tree-level) have been evaluated using the 
FormCalc implementation [28] . 

We will show the results for one representative numerical example. The parameters are 
chosen according to the scenario S shown in Tab. [T|, but with one of the parameters varied. 
The scenarios are defined such that all decay modes are open simultaneously to permit an 
analysis of all channels, i.e. not picking specific parameters for each decay. We will start with 
a variation of rrig, and show later the results for varying ipg. The scenarios are in agreement 
with the MSSM Higgs boson searches at LEP [3ll[35]. The choice of relatively large gluino 
and squark masses also avoids all LHC bounds [36]. The slepton/lepton sector of the cMSSM 
does not enter into our calculation, hence the parameters are not specified. Furthermore, 
also the following exclusion limits [I9] hold in our scenario: 

m^o > 46 GeV, m^o > 62 GeV, m^o > 100 GeV, m^o > 116 GeV, m-± > 94 GeV . (47) 

The scalar quark masses of the parameter set S are shown in Tab. El The values of rrig 
allow copious production of the gluino at the LHC, once ^/s = 14 TeV is reached. 

The numerical results we will show in the next subsections are of course dependent on 
choice of the SUSY parameters. Nevertheless, they give an idea of the relevance of the full 
one-loop corrections. As an example, the various decay widths r(^ — )■ q) are all of similar 
size, contributing similarly to Ftot. The same holds for the various branching ratios. If the 
(artificial) degeneracy of the soft SUSY-breaking parameters in the squark sector were lifted, 
these results could change significantly. However, the size of the loop corrections shown in 
our numerical example still gives a good indication about the size of the expected corrections. 
In the special case of nig < rriq^ (for all q and i = 1,2) the loop-induced decays g Xk9 
{k = 1,2,3,4) as well as three-body decays (which are not investigated here) could become 
dominant. 



11 



Seen. 


tan /3 


Mh± 










Ak 


Ml 


Ms 


Ms 


S 


20 


200 


700 


800 


200 


1000 


800 


200 


300 


1200 



Table 1: MSSM parameters for the initial numerical investigation; all masses are in GeV. 
For the q sector the shifts in M^^ ^ as defined in Eqs. (49) and (50) in Ref. [TB] are taken into 
account. The values for At (= Ac = Au) and A^ (= Ag = Aa) are chosen such that charge- 
or color-breaking minima are avoided [37] . 

















g-type 


u 


c 


t 


d 


s 


b 




697.919 


697.913 


637.268 


702.402 


702.402 


705.879 




799.232 


799.240 


882.565 


800.384 


800.384 


800.848 



Table 2: The scalar quark masses in S for the numerical investigation; all masses are in GeV 
and rounded to one MeV. 



4.2 Full one-loop results for varying rrig 

The results shown in this and the following subsections consist of "tree" , which denotes the 
tree-level value and of "full" , which is the decay width including all one-loop corrections as 
described in Sect. |3l Additionally in this section we also investigate the accuracy of the pure 
SUSY QCD corrections for the decay widths as evaluated for real parameters in Ref. [5]|ll 
The corresponding curves are labelled "SQCD". We start the numerical analysis with g decay 
widths evaluated as a function of rrig starting at rrig = 700 GeV up to rrig = 2 TeV, which 
roughly coincides with the reach of the LHC for high-luminosity running. The upper panels 
contain the results for the absolute value of the various decay widths, r(^ — j- xy) (left) and 
the relative correction from the full one-loop contributions (right), where we compare "tree", 
"full" and "SQCD". The lower panels show the same results (but leaving out "SQCD") for 
BR(^^xy). 

Since all parameters are chosen real no difference between the two decay modes g ^ q\q 
and g ^ (jiq arises, neither in the decay widths, nor in the branching ratios (where of course 
both channels are taken into account in the total decay width). Consequently, we only show 
results for the channel "^ — ?■ g" , and the results for BR(^ ~^ qI Q + 9 ^ Qi q) ^-re simply 
the BR(^ —J- (ji q) multiplied by two. In Sect. 14.31 both channels will be shown separately. 
Furthermore, due to the absense of complex parameters no contributions from absorptive 
parts of self-energy type corrections on external legs can contribute. Again, this will be 
different in Sect. 14.31 

In our numerical scenario S we have chosen a small splitting between the left- and the 
right-handed soft SUSY-breaking parameter within all squark flavors, see Tab. [H If the 
effects of mass differences of the quark and squarks in the final state were neglected one 

^ We have checked that we are in good agreement with Ref. [S] using their input parameters, where a 
small difference remains due to the different renormalization schemes. 
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would expect Sj=i 2BR(^ — )■ g^g) ^ 1/12 (or Sj=i^2 BR(^ Q + 9 ^ QiQ) ~ 1/6 if both 
combinations are taken into account). Deviations from this number are expected due to 
kinematical effects, and differences in the loop corrections (including possible effects related 
to the scalar quark mixing). For large values of rrig where mass effects should be small it is 
expected to reach the value of ~ 1/12. 

We start with the decay g tit. The results for this channel are shown in Fig. [3] as 
a function of nig. The tiny dip in S at nig ^ 1054.6 GeV is the threshold nig = m^^ + mt 
of the self energy S^(m?) in the renormalization constants 6Zg and 6M3. One can see that 
the size of the corrections of the decay width is especially large very close to the production 
thresholdjfl Away from the production threshold relative corrections of ~ —6% are found. 
The difference between the full and the SQCD corrections is roughly 5%; at low masses 
the EW corrections are about one third of the full corrections, while at large nig they are 
dominating over the SQCD corrections. The branching ratio in our numerical scenario S 
(where we have scalar quarks very close in mass) can reach more than 12% close to threshold 
and stays above 5% for large nig. The correction to the branching ratio is then found to be 
~ -4%. 

Next, in Fig. |4] we show the decay r(^ — > t2t). At low nig the decay is kinematically 
forbidden and reaches ~ 11 GeV at nig = 2 TeV, with small positive corrections to the decay 
width at low nig and very small negative corrections at large nig. The behavior of the SQCD 
corrections alone is quite different for this decay. They are large and negative, ~ —14% 
close to threshold and turn positive up to +2% for large nig. The SQCD contributions alone 
constitute a very weak approximation to the full result in this channel. The BR reaches 
nearly 3% at large nig, so that BR(^ — )■ ti t) + BR(^ —7-^2^)^ 8.5% ~ 1/12 is reached. The 
relative corrections reach from 10% where the width is small going down below 1% where 
the width and BR are largest. 

Now we turn to the decays to b/b. The results for the decay g ^ bib are presented in 
Fig. [5l Within S the first dip at nig ^ 800.51 GeV stems from the threshold nig = +Tnc 
in the self energy Eg(m?) entering the renormalization constants 6Zg and The second 

dip at nig ~ 1054.6 GeV comes from the threshold rrig = m^^ + nif. As expected the decay 
width rises from zero at threshold to ~ 15.5 GeV at nig = 2 TeV, including the full one- 
loop corrections. The first threshold leads to very large corrections of nearly —14%, which 
then decreases to ~ —5% at larger nig values. The SQCD corrections approximate the full 
result quite well up to nig = 1200 GeV, but then tend to zero. Again for large nig the full 
corrections are dominated by the EW contributions. The first threshold is also clearly visible 
in the BR, where values around 9% can be reached, going down to ~ 4% at nig = 2 TeV. 
The relative corrections to the BR vary between —4% and +0.5% around the first threshold 
and go to —3% at high nig values. 

The results for ^ — t- 62 ^ are shown in Fig. [6l Again the dip at nig 1054.6 GeV comes 
from the threshold nig = m^^ +nit. The threshold nig = ni^^ +nic and the corresponding large 
effects observed for g bib are absent here due to the fact that = 800.848 GeV. Apart 
from this we find for the decay widths the expected values with r(^ — )■ 62 &) ~ 15 GeV at 

^ It should be noted that a calculation very close to threshold requires the inclusion of additional (non- 
relativistic) contributions, which is beyond the scope of this paper. Consequently, very close to threshold 
our calculation (at the tree- or loop-level) does not provide a very accurate description of the decay width. 
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rrig = 2 TeV (with the SQCD corrections giving a very good approximation to the full result) 
and a corresponding BR of ~ 4%. Again the sum of the two BR's reaches ~ 8.5% ~ 1/12. 
The relative corrections can reach nearly —9% at nig ^ 900 GeV for the decay width and 
values between +9% and +2% for small and large BR's, respectively. 

Next we analyze the effects in the decays involving second generation (s) quarks. We 
start with the decay ^ — ?■ Ci c shown in Fig. [71 Two dips are visible: within S the first 
dip appears at rrig ^ 809.3 GeV, due to the threshold rrig = + mf. The second dip at 
rrtg ~ 1054.6 GeV comes from the threshold rrig = m^^ + rrit. The decay width behaves "as 
expected", rising up to ~ 16 GeV at rrig = 2 TeV. The relative corrections are again largest 
around the first threshold, reaching nearly —10%, flattening out to ~ —3% at larger rrig 
values. The SQCD corrections alone yield larger corrections by ~ —5%, going to zero for 
large rUg, where the EW contributions are dominating. The branching ratio is very large 
at threshold due to the relative smallness of ruc^ + rUc in S, exceeding 20%. At large rUg 
the "expected" ~ 4% are reached. The relative correction to the BR exhibits a maximum 
around rrig = 800 GeV due to rrig = m^^ +mt, see above. At large rrig the corrections become 
very small, but potentially reach larger negative values at rrig > 2 TeV. 

The results for g ^ are shown in Fig. [HI Again, the dip at rrig ^ 1054.6 GeV stems 
from the threshold rUg = rrii^ + rrit. Decay width and branching ratio show the expected 
behavior. The BR reaches ~ 4% at large rrig, so that the sum of the four "charm BR's" goes 
to ~ 1/6. The corrections to the width are maximal with ~ —7.5% around rrtg = 900 GeV 
and rise to small positive values for large rrig, with the SQCD corrections deviating from this 
by up to —1% at most. Concerning the branching ratio the effects are largest where the BR 
is small and reach the level of ~ 3% at rUg = 2 TeV. 

The decays involving scalar strange quarks, ^ — )■ Sj s (i = 1, 2), are shown in Figs. [9l [lOl 
Again several dips are visible, the first dip at rrig 809.3 GeV comes from the threshold 
rUg = rrii^ + rrit. Once more the second dip at rrig ^ 1054.6 GeV is the threshold rrig = 
mi^ + rrit. The results are very similar to the ones for the decays g CiC, where differences 
can only be observed close to threshold due to slightly different masses of the scalar quarks. 
The most prominent example in this respect is the BR(^ — j- si s) that reaches "only" about 
10% for rUg = 700 . . . 800 GeV. Otherwise the "expected" behavior is found at the tree and 
at the loop level, including also the size of the pure SQCD corrections. 

The results for decays involving first generation (s) quarks are very similar to the ones for 
the second generation, again differences appear only due to small deviations in the scalar 
quark masses, see Tab. [2 Otherwise again the "expected" behavior is found with loop 
corrections yielding maximum values of —10% for the decay widths, again with the same 
size of the SQCD corrections, and +5% for the branching ratios. It should be kept in 
mind that the results found for the branching ratios are highly model dependent, and much 
larger /smaller values can be found for different kinematical situations. 
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Figure 3: r(^ -^tit). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 4: r(^ — j- t2t). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 5: r(^ — bib). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 6: r(^ — 62^)- Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 7: T{g ^ Cic). Tree-level ("tree") and full one- loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 8: r(^ — >■ C2c). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 



20 





Figure 9: r(^ — )■ Si s). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 



21 




BR 5BR/BR 




Figure 10: r(^ — )■ §2 s). Tree- level ("tree") and full one- loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 11: r(^ Uiu). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. 
The upper left plot shows the decay width, the upper right plot shows the relative size of 
the corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot 
shows the BR, the lower right plot shows the relative size of the BR. 
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Figure 12: r(^ U2u). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. 
The upper left plot shows the decay width, the upper right plot shows the relative size of 
the corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot 
shows the BR, the lower right plot shows the relative size of the BR. 



24 





Figure 13: r{g ^ did). Tree- level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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Figure 14: r(^ — )■ ^2 (i). Tree- level ("tree") and full one-loop ("full") corrected decay widths 
are shown with the parameters chosen according to S (see Tab. [1]), with rrig varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. Also shown are the pure SQCD corrections ("SQCD"). The lower left plot shows 
the BR, the lower right plot shows the relative size of the BR. 
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4.3 Full one-loop results for varying cpg 

In this subsection we analyze the various decay widths and branching ratios as a function of 
i^g. The other parameters are chosen according to Tab. [H 

When performing an analysis involving complex parameters it should be noted that the 
results for physical observables are affected only by certain combinations of the complex 
phases of the parameters /i, the trilinear couplings At, A^, . . . , and the gaugino mass param- 
eters Ml, M2, M3 [25l[38]. It is possible, for instance, to rotate the phase ipM2 away, and we 
remain with ipM^ and v'Ms = as independent physical phases. Experimental constraints 
on the (combinations of) complex phases arise in particular from their contributions to elec- 
tric dipole moments of heavy quarks |39] , of the electron and the neutron (see Refs. [^BT] 
and references therein), and of deuteron [42], Refs. [40ll43ll44] for reviews. The phase of 
the gluino mass parameter remains relatively weakly constrained. Since our decaying par- 
ticle is the gluino we focus here on the dependence on ipg and set all other phases to zero. 
Analyses of the full one-loop effects of SUSY decays on ipA^, ^At and ifu^ can be found in 
Refs. [T5l[T6l[271, respectively. 

Since now the complex phase of M3 can appear in the couplings, two effects arise: first a 
difference between r(^ — )■ q\q) and r(^ — )■ g) is observed, and consequently we now show 
the results separately for these two decay channels (but stick to the notation — )■ gj g" when 
referring in general to a plot). These differences arise naturally from the combination of an 
imaginary (absorptive) part of a decay diagram and a complex coupling. As required, we 
always find r(^ — )■ q\ q)\ip- = T{g — )■ g^ q)\~ipg- A second source of this type of contributions 
comes from absorptive parts of self-energy type corrections on external legs (called "absorp- 
tive self-energy contributions" from now on), and they also have been included according 
to the formulas given in Sect. 12. 3[ For (fg = 0, vr, 2-7r these "absorptive effects" vanish (by 
construction) . 

As before we show the decays in Fig. [15] - [26] The arrangement of the panels is the 
same as in the previous subsection (but leaving out the pure SQCD corrections). According 
to Tab. [T] we choose rrig = 1200 GeV, which is substantially larger than the squark masses. 
Accordingly we expect again that BR(^ — gi g) + BR(^ — g2 g) ~ 1/6, with possible 
deviations in the decays involving scalar tops, where the mass variation is maximal. 

The results for g ^ tit are shown in Fig. [TS] where the results are given as a function 
of ipg. One can see that both the decay widths as well as the size of the corrections to it 
vary substantially with ipg. While the widths ranges between ~ 8 GeV and ~ 4 GeV, the 
corrections vary between —8.5% and —2.5% in our numerical scenario S, where also a small 
difference between g i\t and g iii can be observed. The BR's vary correspondingly 
from ~ 8.0% to ~ 4%, where the size of the loop effects in S goes from —4% to ~ -|-2%, again 
with a small visible difference between g ^ i[t and g ^ iii. The second channel involving 
scalar tops, ^ — ^2 ^, is shown in Fig. [TH] In agreement with Fig. [1] the decay width is 
substantially smaller than for g tit, largely due to the reduced phase space, varying 
between ~ 1 GeV and ~ 3.5 GeV. The difference between r(^ — ilt) and r(^ — 
is clearly visible. The loop corrections range from -|-4% to —10%, with a clear symmetry 
between g ^ i^t and ^ — )■ ^2 ^- Correspondingly, the BR(^ — )■ ^2 t) reach values between 
~ 1% and ~ 3.5% with a variation of -|-8% to —6%. The sum of the four BR's is close to but 
always a bit smaller than ~ 1/6 due to the non- negligible effects of the scalar top masses 
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and mixings. 

For all other decay modes, i.e. g QiQ with q ^ the mass and mixing effects are 
substantially smaller|j Within S we find T{g ^ qi q, q ^ t) (shown in Figs. [T71 Uni [211 [231 
[25]) at the level of ~ 5.5 GeV, where small variations depend on the squark flavor. Also the 
required symmetry of r(^ QiQyQ 7^ t)\fg ~ ^id ~^ QiQyQ t)\-fg is clearly visible, and 
similarly for the BR's and the size of the loop corrections. While the size of the one-loop 
corrections vary between —1% and —6% in the case of ^ — ?■ 61 b, they are found to be within 
~ —3% and ~ —4.5% for the other flavors. The corresponding BR's are found at the level 
of ~ 5% with less than ±0.25% of variation. The relative corrections in S are at the few 
per-cent level at most. 

The results for the decay modes to the heavier scalar quarks, g ^ q2(l,(l ^ t, shown 
in Figs, [m [20I [221 [Ml [261 cire also very similar, again with a visible difference between 
g ^ qIq and g ^ q2q,q t. The widths are ~ 4 GeV with a small phase dependence. 
The size of the one-loop corrections vary around —5% for ^ — >■ 62 ^, ~ —4.35% for g — )■ 
= c, u, and ~ —4.85% for g ^ q2q,q = s,d. The BR's for all flavors are found at 
~ 3.5% with loop corrections at the per-cent level. As "expected", for all flavors we find 
BR(^ ^ gi g) + BR(^ ^ g) ~ 1/6. 

4.4 The one- loop decays g ^ X^g 

In Figs. [27|-[30] we present the variation of r(^ — x° g), {k = 1, 2, 3, 4) as a function of rrig 
and ifg in S, see Tab. [H The structure at rrig ^ 1054.6 in Fig. [271 and Fig. [28] is the vertex 
production threshold rrig = m^^ + m^. In Fig. [221 and Fig. [SU] the dip at rrig ^ 699.2 GeV 
stems from the vertex production threshold rrig = rric^ + m^, the peak at rUg ^ 809.3 GeV is 
the vertex production threshold rUg = + m^. 

With our choice of parameters, see Tab. [H we find the neutralino masses at ~ 151, 206, 
218, 338 GeV, i.e. relatively similar with respect to the larger gluino mass. The size of the 
decay widths is very small for all four neutralinos with respect to the hadronic final states 
discussed in the previous subsection. The largest values are reached for g — )■ y^g, going up to 
~ 0.007 GeV. The dependence on rrig is very similar, as expected from the large gap between 
neutralino and gluino mass. We see a threshold on-set, going up to the maximum values 
around 800 GeV. Due to the smallness of the decay widths the corresponding branching 
ratios are negligible in However, in a situation where the gluino is lighter than all scalar 
quarks these loop-induced two-body decay modes can become relevant and may constitute 
a significant source for x° production at the LHC (which could be used to measure the 
properties of the SUSY CDM candidate). In such a situation, of course, also three particle 
final states will be relevant, which are, however, beyond the scope of this paper. 



^ Especially the tree- level results depends strongly on the scalar quark mass splitting rriqXq (sec Eqs. (|4]), 
dni) via C/q^je^^s/^ and is therefore nearly constant with respect to ipg for q ^ t. 

They arc not included in the calculation of the total decay widths shown in Sect. 14.51 
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Figure 15: r(^ — )■ tit). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 16: r(^ — )■ t2t). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 17: r(^ — j- bi h). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 



31 



r/GeV 












4.1 






4.4% 


1 


1 1 1 


1 


1 1 


4.05 






-4.6% 
-4.8% 




\ 






4 




62 tree 

bib, full ■ 


-5% 


/ 


\ / 






3.95 




k k M 


-5.2% 


/ 






/)t Jj , \ 

hi / 


3.9 






-5.4% 
-5.6% 








/ 


3.85 






-5.8% 










3.8 




1 1 1 1 1 


6% 


1 


1 1 


1 





0° 45° 90° 135° 180° 225° 270° 315° 360° 



0° 45° 90° 135° 180° 225° 270° 315° 360° 




Figure 18: r(^ — > &2^)- Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 19: r(^ — Cic). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 20: r(^ — C2c). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 21: r(^ — Si s). Tree- level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 22: r(^ — S2s). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 23: r(^ Uiu). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 24: r(^ — U2u). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 25: r(^ — )■ did). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 26: r(^ — )■ d2d). Tree-level ("tree") and full one-loop ("full") corrected decay widths 
(including absorptive self-energy contributions) are shown. The parameters are chosen ac- 
cording to S (see Tab. [1]), with ipg varied. The upper left plot shows the decay width, the 
upper right plot shows the relative size of the corrections. The lower left plot shows the BR, 
the lower right plot shows the relative size of the BR. 
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Figure 27: r(^ — Xifi')- One-loop ("loop") corrected decay widths are shown with the 
parameters chosen according to S (see Tab. [1]), with rrig (left plot) and ^pg (right plot) 
varied. 
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Figure 28: r(^ — x^g). One-loop ("loop") corrected decay widths are shown with the 
parameters chosen according to S (see Tab. [1]), with rrig (left plot) and tpg (right plot) 
varied. 
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Figure 29: r(^ — Xsfi')- One-loop ("loop") corrected decay widths are shown with the 
parameters chosen according to S (see Tab. [1]), with rrig (left plot) and ^pg (right plot) 
varied. 
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Figure 30: r(^ — x^q)- One-loop ("loop") corrected decay widths are shown with the 
parameters chosen according to S (see Tab. [1]), with rrig (left plot) and tpg (right plot) 
varied. 
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4.5 The total decay width 

Finally we show the results for the total decay width of the gluino, see Eqs. ( HTj) . In Fig. |31] 
the upper panels show the absolute and relative variation with rrig. We observe a nearly linear 
rise from the threshold region on up to rrig = 2 TeV, where Ftot ~ 370 GeV is reached. The 
relative size of the corrections are around +15% close to the production threshold, reaches 
—10% at nig = 800 GeV and levels out at ~ —2.5% for large gluino masses. Depending on 
the choice of SUSY parameters and the location of (production) thresholds the size of the 
one-loop corrections are non-negligible at the LHC. We also included here the pure SQCD 
corrections, which tend to overestimate the full result by ~ —5% at low rrig and go to zero 
at large rrig, where the EW corrections are dominating. 

The lower panels of Fig. [3T] show the result as a function of (pg. For the nominal value 
of rrig = 1200 GeV we find the loop-corrected total decay width between ~ 113 GeV for 
(Pg = 0° and ~ 110 GeV for (fg = 180°. The relative size of the corrections vary between 
—4.8% and —4.2%, possibly still relevant for LHC measurements. 

5 Conclusions 

We have evaluated all two-body decay widths of the gluino in the Minimal Supersymmetric 
Standard Model with complex parameters (cMSSM). The decay modes are given in Eqs. ([1]) 
- ([3]). After the (pair) production of gluinos they are expected to decay via cascades to 
lighter SUSY particles and quarks and/or leptons. Consequently, the decay of a gluino to a 
scalar quark and a quark will one of the most relevant decays in such a cascade. In order to 
determine the masses, couplings and parameters of SUSY particles precisely in such cascades, 
higher-order corrections are necessary to keep the theory error on an acceptable level. 

Our evaluation is based on a full one-loop calculation of all decay channels, also including 
soft and hard QED and QCD radiation, necessary to derive a reliable prediction of any two- 
body branching ratio. We also include the purely loop induced decays ^ — Xfcfi' (^ = 
1, 2, 3, 4). These modes as well as possible three-body decay modes can become sizable only 
if all the two-body channels are kinematically (nearly) closed. We have not investigated 
these three-body decay modes, which are, beyond the scope of this paper. 

We first reviewed the relevant sectors of the cMSSM, including their renormalization 
(where extensive explanations can be found in Refs. [ini[ISl[27] ) . We have discussed the 
calculation of the one-loop diagrams, the treatment of UV- and IR-divergences (where the 
latter are canceled by the inclusion of real soft QCD and QED radiation), and in addition 
the hard QCD and QED bremsstrahlung. 

For the numerical analysis we have chosen one parameter set that allows simultaneously 
all two-body decay modes and respect the current experimental bounds on Higgs boson 
and SUSY searches. In a first step we have shown the results as a function of rrig. For 
the largest value included in our analysis, rrig = 2 TeV we find decay widths of ~ 15 GeV, 
and BR(^ — >■ qiq) + BR{g — >■ q2q) ~ 1/6. The largest deviations were found for the 
decays involving the third generation (s) quarks, where effects of the squark mixing and 
corresponding changes in the masses and couplings are maximal. The size of the one-loop 
corrections is found to be largest close to the production thresholds and reaches several 
per-cent for large rrig. We have also investigated the quality of the approximation of the 
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Figure 31: Ttot- Tree-level ("tree"), full one-loop ("full") and pure SQCD ("SQCD") cor- 
rected decay widths are shown with the parameters chosen according to S (see Tab. [1]). The 
upper left plot shows the total decay width, the upper right plot shows the relative size of 
the total corrections, with rrig varied. The lower plots show the same but with ifg varied 
(including absorptive self-energy contributions). 
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pure SQCD corrections. While for some channels, especially involving lighter quarks, the 
SQCD contributions yield a good approximation, in the channels involving scalar tops and 
bottoms the EW corrections, which were evaluated here for the first time, are crucial to yield 
a reliable result. 

In a second step we have analyzed the dependence of the decay widths and branching 
ratios on the phase of the gluino mass parameter, ipg, for = 1200 GeV fixed. Again, the 
results are to a large extent flavor independent, except for decays involving (scalar) tops. 
In the latter case a variation of the BR of up to 14% was found, and the size of the loop 
corrections can yield ~ 10%. For the other flavors we flnd as before BR(^ — ^ Q') + BR(^ — 
q2q) ~ 1/6, and a smaller variation with ipg. 

The size of the widths of the purely loop-induced decays is substantially smaller, staying 
below 0.01 GeV for our parameters. This could change if the gluino were lighter than all 
scalar quarks. In that case also three-body decays would be relevant, which, however, were 
not investigated here. 

Furthermore, it has to be kept in mind that we have evaluated the full one-loop corrections 
for one set of parameters only, which was chosen to have many decay modes open, and not to 
emphasize the loop corrections. Consequently, it can easily be imagined that substantially 
larger or smaller corrections can be found for other choices of SUSY parameters. The full one- 
loop corrections should (eventually) be taken into account in precision analyses of cascade 
decays at the LHC Nevertheless, the size of the loop corrections found give an idea of the 
relevance of the full one-loop corrections. 

Following our analysis it is evident that the full one-loop corrections are mandatory for 
a precise prediction of the various branching ratios. The results for the gluino decays will 
be implemented into the Fortran code FeynHiggs. 
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